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We demonstrate that the gauged BPS baby Skyrme model with a double vacuum poten¬ 
tial allows for phase transitions from a non-solitonic to a solitonic phase, where the latter 
corresponds to a ferromagnetic liquid. Such a transition can be generated by increasing the 
external pressure P or by turning on an external magnetic held H. As a consequence, the 
topological phase where gauged BPS baby skyrmions exist, is a higher density phase. For 
smaller densities, obtained for smaller values of P and H, a phase without solitons is reached. 

We hnd the critical line in the P, H parameter space. Furthermore, in the soliton phase, we 
hnd the equation of state for the baby skyrmion matter V = V{P, H) at zero temperature, 
where V is the ’’volume”, i.e., area of the solitons. 

PACS numbers: 


I. INTRODUCTION 

The baby Skyrme model [T] (see also m-m) is a lower dimensional version of the 3-1-1 dimen¬ 
sional Skyrme model [TTj, that is, an effective theory of the low energy, nonperturbative sector of 
Quantum Chromodynamics. It plays an important role as a toy model for the higher dimensional 
theory. In fact, many aspects of the original Skyrme proposal were understood qualitatively by 
the investigation of its 2 -|- 1 dimensional counterpart. Recently, it has been understood that the 
correct Skyrme model should be a near BPS theory m, m where the BPS part gives the main 
contribution to some static observables m (see also [18], |19ji. A concrete near BPS Skyrme 
action has been proposed in US). Indeed, already its BPS part leads to very accurate results not 
only for the binding energies of the atomic nuclei [T5| and some thermodynamical properties |20| 
of the nuclear matter but also for the physics of the neutron stars [2T] (maximal masses and radii, 
equation of state). There are, therefore, strong arguments for the assumption that this near BPS 
Skyrme model, providing an unified description of the nuclear matter from baryons with the baryon 
charge B = 1 and nuclei B ~ 10^ to neutron stars B ~ 10^^, is a good candidate for the correct 
low-energy effective model of QCD. 

Qualitatively, the BPS Skyrme model consists of two terms: the topological (baryon) current 
squared and a potential. Obviously, there is a baby version of the BPS Skyrme model [22], [23] . 
|24| . |25j (for some near BPS extension see [26]). Both theories have many features in common: 1) 
the existence of the linear energy - topological charge bound which is saturated by infinitely many 
solutions with arbitrary topological charge; 2) the invariance of the static energy integral under the 
SDiff symmetries and the resulting close relation to an incompressible liquid as the corresponding 
state of matter; 3) the existence of a standard thermodynamics with equivalence between thermo¬ 
dynamical and field theoretical pressure. 
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Recently, the baby BPS Skyrme model has been used to investigate magnetic properties of 
skyrmions after minimal coupling to the Maxwell electrodynamics m, [28] . Interestingly, the 
gauged baby BPS Skyrme model also allows for a = 2 supersymmetric extension, see, e.g., [29], 
m,m- Undoubtedly, the interaction with an abelian gauge field is very important for properties 
of baryons and atomic nuclei. It is a well-known result by Witten that, although the precise form 
of the low energy effective Skyrme-type theory is still unknown, its coupling to the electromagnetic 
C/(l) field is completely fixed |3T]- Hence, in principle, one can study electromagnetic properties 
of skyrmions once a particular model is chosen. In practice, however, this is a rather complicated 
task even within the framework of the BPS Skyrme model, and not many results are available (for 
the standard Skyrme model see [32]-[33| and [3l|). Again, the BPS baby Skyrme model serves here 
as a perfect laboratory (for gauged solitons in a related Skyrme-Faddeev-Niemi model see [35] i. 
The aim of the present paper is to further investigate the gauged BPS baby Skyrme model EZj, 
[28] . Specifically, we would like to analyze the issue of the (non) existence of planar skyrmions in the 
model with double vacuum potentials. It is one of the most surprising findings in our former work 
that, in contrast to one vacuum potentials, the gauged BPS baby Skyrme model does not support 
solitons for potentials with two vacua. This unexpected observation (there is no such difference 
for the non-gauged case, see for example recent works for the Skyrme model |36|) was made for 
BPS configurations, i.e., with zero external pressure and no external magnetic field. Here we will 
show that solitons can exist if these two external parameters are larger than some minimal eritical 
values. Hence a phase transition between a solitonic phase and a non-solitonic phase is observed. 


II. THE GAUGED BPS BABY SKYRME MODEL 


The gauged BPS baby Skyrme model is given by the following Lagrangian density 

where the covariant derivative reads m, 138] 

+ A^n X (H.2) 

As usual, without loss of generality, we assume that the constant vector n = (0,0,1) and the 
potential C/ is a function of the third component of the unit vector field (p = (cpi, <p2, (ps)- The 
gauge field is the Maxwell field with the electromagnetic coupling constant g. Here we list the 
main properties of this model. 


i) BPS bound 

The static energy (or regularized energy for a non-zero constant external magnetic field H) 

Eo 


Ereg — 


d X 


X D2$y + 2fi^U + 4(^ - 


(H.3) 


is bounded from below by the topological charge k 


Ereg > 47r|A:|UoA^ < W' >S2 . 


(11.4) 
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Here B is the magnetic field (pseudoscalar in 2+1 dimensions), W is the so-called superpo¬ 
tential (depending on the target space variable (ps) dehned as 

- 2A^WH = 2fi^U, (II.5) 

and < W' >52 is its average value over the target space manifold. The prime denotes 
differentiation with respect to 

a) BPS equations 

The bound is saturated by solutions of the BPS equations 

Q = W' (II. 6 ) 

B = -g^A^W + H (II.7) 


where 


Q 


q + eijAidj{n ■ q = cf ■ di$ x d 2 (j). 


(II. 8 ) 


and q is the topological charge density. 

Hi) No solitons for two-vacuum potentials 

The existence of such solutions is a rather subtle issue as the superpotential W has to be 
dehned over the whole interval 4>3 G [—1,1]. This requirement is very restrictive and, for 
example, cannot he satisfied in the case of two vacuum potentials and H = 0. The fact that 
the BPS baby skyrmions (two vacuum potential) can be destabilized by the electromagnetic 
interaction was quite surprising since no similar effect was known in the non-gauged case. 


iv) Magnetization 

The thermodynamic magnetization M, dehned as minus the change of the thermodynamic 
energy of a sample (in our case, the skyrmion) under a variation of the external magnetic 
held 


M = - 


dE, 


reg 


dH 


(II.9) 


exactly agrees with the usual dehnition of magnetization as the difference between full and 
external magnetic hux in the sample 



jiB- H)d^x. 


(II.IO) 


Furthermore, the matter described by the gauged BPS Skyrme model behaves as a rather 
nonlinear ferromagnetic medium where the solitons remain magnetized even when the ex¬ 
ternal magnetic held vanishes, i.e., they possess a permanent magnetization. The external 
magnetic held H acts in two ways: it change the size of the baby skyrmions (it squeezes 
solitons if has the same sign like the permanent magnetization M of the skyrmion, while it 
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enlarges the skyrmion if H and M have opposite signs) and interferes with the magnetic 
field inside the solitons: for an external magnetic field H which is sufficiently weak and 
oppositely oriented to M, the phenomenon of magnetic flux inversion occurs. That is, the 
total magnetic field B flips sign in a shell region near the boundary of the skyrmion, as it has 
to take the value B = H at the boundary. On the other hand, it preserves its original sign 
resulting from the permanent magnetization in the interior (core region) of the skyrmion. 

v) Non-BPS (external pressure) solutions 

Solutions with a non-zero field theoretical pressure {P = —Tu, where is the energy- 
momentum tensor with the constant external magnetic part subtracted) are described by 
the first order BPS equations but with a modified superpotential equation 

- 2\^WH = 2p,^U + 2P. (II.ll) 


The pressure introduced in this way obeys the usual thermodynamical relation 

p = (ii. 

where V is the geometrical (compacton) ’’volume” (area) and therefore agrees with the 
thermodynamical pressure. 


III. THE NEW BABY POTENTIAL 


In order to understand how the external magnetic field and/or pressure may influence the 
existence of planar solitons in the case of two vacuum potentials we have to analyze the necessary 
condition for the existence of BPS gauged baby skyrmions, that is, the existence of a global solution 
to the superpotential equation 

A^W'^ + g^A^W^ - 2A^WH = 2g^U + 2P. (III.l) 


Concretely, we shall consider the so-called new baby potential 

{ 111 . 2 ) 

The next step (if the superpotential exists) is to solve the corresponding equations of motion. Here 
we assume the standard axially symmetric static ansatz 


(f{r, 4>) 


( sin f (r) cos ncj) \ 
sin f{r) sin ncj) , 
\ cos /(r) / 


Aq = Ar = 0, A^ = na{r) 


(III.3) 


which leads to an identically vanishing electric field and to the magnetic field B = Moreover, 

we introduce a new base space variable y = r^/2 and target space variable h via 


(j)^ = cosf = l-2h /i =-(1 - cos/). 


hy = - sin ffy. 


(III.4) 
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Then the first order field equations (BPS equations) read 

2nhy{l + a) = --Wh (III.5) 

nay = -g^X^W + H (III.6) 

where the following boundary condition 

h(0) = 1, h{yp) = 0, (III.7) 

a(0) = 0, ay{yp) = ^ (III.8) 

guarantee the nontrivial topology (n is just the topological charge) and finiteness of the (regular¬ 
ized) energy 

Ereg = 27r J (ly ^2A^n^(l a)‘^hy + y^h{l - h) + • (III.9) 

Here, yp is the size of the compacton for a non-zero value of the pressure P. Since we deal with 

compact solutions one can define the geometrical volume (area) 

V = 2TTyp (III.10) 

which, as we remarked before, together with the pressure fulfills the usual thermodynamical rela¬ 
tion. The superpotential equation in terms of the new variable h reads 

\2 

—W^ + c/^A^IT^ - 2X^WH = 2y^U + 2P. (III.ll) 

Further, we always assume that the Skyrme field approaches the vacuum at /i = 0 at the compacton 
radius rp = \/2yp. This, together with the boundary conditions and BPS equations for h and a, 
implies that the superpotential W{h) must obey the boundary conditions 

tT(0) = 0, IT/i(0) = (III.12) 

A 

For nonzero pressure this implies that the first derivative hy jumps at the compacton boundary. 
This is as it should be, because if we impose a nonzero external pressure at the soliton surface 
(compacton boundary), then the energy density must jump there. 

In numerical computations we assumed A = 1 and y? = 0.1, while the electromagnetic coupling 
constant was kept as a free parameter. Furthermore, a particular value of the topological charge 
n can be eliminated from the field equation by the coordinate redefinition y = y/n rendering the 
numerics n independent (with a linear dependence of the energy and volume on n - as expected 
for a BPS model). Therefore we present all plots for n = 1. 
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FIG. 1: The position of the singularity in the superpotential W ior H = Q and P = 0, as a function of the 
coupling constant g. 


A. iJ = 0,P = 0 

Let us begin with the case where LT = 0, P = 0. It is straightforward to see that we arrive at 
some difficulties. Indeed, the superpotential equation 


^Wl + g^\V = 2^?h{l-h) 


(III.13) 


enforces boundary conditions at both ends of the unit interval. In addition to W{0) = 0, we 
must impose W{1) = 0. However, for a first order equation it is in general impossible that both 


conditions can be simultaneously obeyed. Indeed, if we numerical solve equation (III. 13) then a 
singularity occurs for h = Hq. (Strictly speaking Whiho) = 0 while Whh{ho) = oo.) It approaches 
the boundary of the interval as the electromagnetic coupling constant goes to 0, see Fig. 1. Then, 
for = 0 we arrive at the non-gauge case for which BPS skyrmions do exist and read 


Hy) = 


sin^ ^ (l — ^ 
2 V yo 

0 


y <yo 
y>yo 


where yo = 


mrX 


(III.14) 


Obviously, the nonexistence of the superpotential shows that gauged solitons do not show up if the 
pressure as well as the external magnetic field vanish (for any value of 5 > 0). It is also clearly 
visible that in the more general situation, with nonzero value for the external parameters P and 
H, the superpotential might exist on the full unit interval for some external parameter values. 

B. P = 0,P>0 


Now let us turn on the pressure, still assuming H = 0, at the boundary of compactons. For the 
limiting case, when g = 0, the exact solutions are known [28] 

1-2/1 , g , ^ u . yyo g 

2^^ + h{l-h) n^/ 2 A VP 


(III.15) 
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FIG. 2: The position of the singularity for the superpotential for H = 0. For ho > 1 solutions do exist. 


leading to the following equation of state 

4^/271 Xn /i 

V = -arctan . 

^ Vp 


(III.16) 


Hence, for P —)• 0 we find the maximal volume Vmax = ■ For higher charge solutions 

the volume can be further increased by taking a collection of smaller charge constituents with an 
empty space (vacuum) in between. Such a phase resembles a gas of BPS skyrmions. On the other 
hand, for V < Vmax we have a sort of liquid phase. Therefore, one can conclude that we observe 
a gas-liquid phase transition at P = 0 (and no electromagnetic coupling). In fact, such a behavior 
is a generic feature of the BPS Skyrme models in any dimensions [20) . 

If the electromagnetic coupling constant is not zero then solutions exist only for P > Pmin- At the 
critical minimal pressure Pmin the system reaches its maximal volume Vmax{P)- For P < Pmin the 
superpotential develops a singularity for ho < 1, see Fig. 2, and solutions cease to exist. Hence, 
at Pmin another phase transition occurs between a topologically trivial phase (no solitons) and a 
skyrmionic matter phase. One can also say that the gauged BPS skyrmions are pressure induced 
objects. In Fig. 3, we plot the equation of state V = V(P) for the skyrmionic phase for several 
values of the electromagnetic coupling constant. The ends of the curves correspond to Pmin and 
Vmin- A detailed plot is presented in Fig. 4, where the line of minimal pressure is visible as a 
function of g. 


C. i7^0,P = 0 

Another possibility to create solitons, i.e., to extend the existence of the superpotential on the 
whole unit interval, is to consider a non-vanishing external magnetic field P > 0 together with 
P = 0, 

+ g^X^W^ = 2g^U + 2X^WH. 


(HI.17) 
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FIG. 3: Volume vs. pressure, for different values of the coupling constant g. 



FIG. 4: The volume as a function of pressure P and coupling constant g. The white region denotes no 
solutions. 


Let us start with the limiting case when g = 0. The corresponding superpotential equation is 

\2 

— Wl = 2p^U + 2X^WH (III.18) 

It is not difficult to show that Wh will never take the value zero on the unit segment (except 
at /i = 0). Indeed, assuming the positive square root we find that Wh > 0 in the vicinity of 
h = 0. Then, the right hand side of the equation is always greater than 0 as 1/ > 0 and W is 
a growing function. A consequence of that is the existence of solitons for any H (for g = 0 and 
P = 0). Again, if 51 > 0, the original {H = 0) singularity moves towards the boundary of the 
unit segment, reaching it for a certain minimal value of the magnetic field Hmin, Fig. 5. Starting 
with this minimal value, the external magnetic field creates a skyrmionic phase. Now, the volume 
increases with the magnetic field. Fig. 6 , as the magnetization inside the solitons is opposite to 
the orientation of the external filed. A more detailed equation of state V = V{H) for different g is 
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FIG. 6: Volume vs. pressure, for different values of the coupling constant g. 


presented in Fig. 7. 

It is a surprising feature of the gauged BPS baby Skyrme model that the limiting case, g = 0, 
can be solved exactly although the superpotential is not known analytically. Let us write the full 
second order field equation for the Skyrme field assuming the non-back reaction approximation 
(equivalent to (7 = 0), i.e., a = f3y where (3 = H/n. Then, 

dy [hy{l -h a)^] = ^^2)^2 ^ ~ ~ (IIL19) 

Now we introduce the new variable 


Then, 


1 

i + /3y' 


(111.20) 


z hzz — 




4n2A2/32 


(1 - 2h) = a{l - 2h) 


(III.21) 
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FIG. 7: The volume as a function of the external magnetic field H and coupling constant g. The white 
region denotes no solutions. 


with the following general solution: for a < | 


x/U^+l -x/U^+1 1 

h{z) = Ciz 2 +C 2 Z 2 + 


for a > I 


, / X /- ■ (- 1 ^ \ /- ( V 8 a - 1 \ 1 

ri[z) = Liyz sin I -^-In z 1 + G 2 v ^ cos I -^-In z I + -. 


The boundary conditions 

h{z = 1) = 1, /i(zo) = 0, hz{zo) = 0 
lead to the following solution (a < 1/8) 


h{z) = 




1 


(III.22) 


(III.23) 


(III.24) 


+ cot(wlnzo)^ sin(wlnz) + -\/zcos(a;lnz), (III.25) 


2 2 \yTosin(a;lnzo) 

where oj = y/8a^^ and the position of the compacton boundary is given by 

sin(u; In zq) + co cos{uj In zq) = —2a;-^/zo. 

Hence, we find the equation of state relating the volume V = 27ryo and the magnetic field 


(III.26) 


sm 


2 V 2i?2A2 




— 1 In ( 1 + 


HV\ 

Inn J 


+ 2 




2H2\2 


— 1 cos 


2 V 2H2\2 




— 1 In ( 1 + 


ffV'). 

27rn J 


= -2 


_1 

2//2A2 C 



HV 

27rn 


(III.27) 
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Hr, 


FIG. 8: Critical line for different values of the coupling constant g. Above the line the topological phase 
exists while below the line we have the non-topological phase. 


for the external magnetic field obeying 


(111.28) 

If the magnetic field is bigger than this value one has to change to the hyperbolic functions resulting 
in 


— sinh 


2 V 2iJ2A2 


In 1 + 


HV 

27rn 


+ 2 a 1 - 




2H^X^ 


cosh 


-W1 - 

2 V 2^2 A2 




In 1 + 


HV 

27rn 


= -2 




1 H 

l+HV_ 
^ ^ 2 ™ 


(III.29) 


These results give us the approximated equation of state V = V{H) for small value of the electro¬ 
magnetic coupling constant. 


D. H^0,P^0 

In the most general situation, one imposes the external pressure as well as the external magnetic 
field. Again, solitons are observed for sufficiently large P and H. In fact, one can plot a critical 
curve in the P, H plane above which the solitonic phase exists, see Fig. 8. The corresponding 
equation of state V = V(H,F), for g = 0.2, is presented in Fig. 9. The dependence of the 
regularized energy, as well as its derivative, on the magnetic held H for several values of P (and 
g = 0.2) are presented in Fig. 10. 

The soliton medium possesses a permanent magnetization, that is, a non-trivial magnetic held 
inside the solitons, even for vanishing external magnetic held. Of course, this makes sense only if 
such solitons do survive in the H = 0 limit, which requires a sufficiently large non-vanishing external 
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FIG. 10: Regularized energy and its derivative as a function of H for different value of the pressure P. Here 
5 = 0.2 


pressure. We plot the magnetization in Fig. 11. It is always negative and grows to a constant 
value in the limit of very large pressure. In Fig. 12 we show the magnetic susceptibility at H = 0. 
There is a qualitative different to the one-vacuum potential case: although the regularized energy 
is a monotonously growing function with H (positive hrst derivative and negative magnetization), 
it is a concave function (with negative second derivative respect to H). This happens for any P. 
Hence, the magnetic susceptibility (at H = 0) is always positive. 


Xo{P) = 


d'^Ereg 

dH^ 


> 0 . 

H=0 


(III.30) 


This behavior does not change the type of the magnetic medium described by the model. It is still 
a ferromagnetic medium as we have a permanent magnetization, and the field wants to join the 
external held smoothly. It simply shows that the BPS Skyrme model with the new baby potential 
responds in a very nonlinear way to the external magnetic held. 
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FIG. 12: Susceptibility at i/ = 0 as a function of the pressure, for different g 

IV. SUMMARY 

In this paper, we have analyzed the gauged BPS baby Skyrme model with a two vacuum 
potential. It has been shown that although there are no baby skyrmions in the ’’unperturbed” case 
(no pressure and vanishing magnetic field at the compacton boundary), such solitons do appear if 
an external pressure and/or a constant external magnetic field are applied. Strictly speaking, baby 
skyrmions exist if these parameters are above their critical values Pmim Hmin- The resulting critical 
curves (in the P, H plane) have been obtained for several values of the electromagnetic coupling 
constant. Such a curve divides the system into two phases, a solitonic and a non-solitonic one, and 
crossing the curve corresponds to a phase transition from the topological to the non-topological 
phase. 

Otherwise, the thermodynamical and magnetic properties of the model in the solitonic phase are 
qualitatively very similar to the previously analyzed one vacuum case. The baby skyrmions form a 
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FIG. 13: Average density for g = 0.1. The white region denotes the non-solitonic phase. 

non-linear ferromagnetic, i.e. permanently magnetized, type of matter (except for the P < Peru id) 
case, where the non-solitonic phase is reached as H —)> 0). Let us note that a topology changing 
phase transition is a known phenomenon in Skyrme models. For example, a transition between 
phases with skyrmions and half-skyrmions was proposed in the context of dense nuclear matter 
j39| (neutron stars). Specifically, the crystal phase formed by skyrmions transforms, as the density 
of the matter increases, into a different crystal phase where the building blocks are half-skyrmions, 
instead. In our case, however, the approach to the phase transition is quite different. Namely, one 
can start in the topological phase and then, assuming a fixed value of the external magnetic field 
H, reduce the pressure until one reaches the non-topological phase. As the pressure decreases, the 
size of the skyrmions grows (pressure always squeezes the solitons in the BPS Skyrme models), 
and the average energy density {paverage = EreglV) decreases. Fig. 13. Hence, the non-skyrmion 
phase occurs for lower densities. In other words, in the model with the new baby potential, the 
BPS baby skyrmions appear only if the medium has a sufficiently high density. 

There are several directions in which the present investigation may be continued. For example, 
one may be interested in how the magnetic and thermodynamic properties of the gauged BPS 
baby skyrmions are modified by the (iso)spin. The isospinning BPS baby skyrmions are known 
|22j . and it would be interesting to look for such solutions also in the gauged version of the model 
(Isospinning baby skyrmions and 3-|-l dim skyrmions have been recently considered in [8] and |40j ). 

Unfortunately, there are a lot of problems and modifications which must be solved before we 
could extrapolate these results to the 3-1-1 dimensional (near) BPS Skyrme model. First of all, the 
planar case is a bit special with respect to the magnetic field. Indeed, here i? is a pseudo-scalar 
instead of a pseudo-vector. Moreover, in contrast to to the planar case, the gauging of the BPS 
Skyrme model will break the SDiff invariance in the static sector. Furthermore, the proper, QCD- 
induced, coupling to the abelian gauge field requires also the Wess-Zumino-Witten term, which 
leads to a non-zero electric field inside the skyrmions. Finally, the inclusion of non-BPS terms, i.e., 
the sigma model term and the Skyrme term, may modify the magnetic properties of the system. 
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In any case, the first step should be a derivation of the energy minimiser in the non-gauge near 
BPS Skyrme model, which is a difficult numerical issue |H] (for some recent numerical result in 
the Skyrme model with the sextic term see [36], [42| 1. 

Moreover, there exists the possibility that the gauged BPS baby Skyrme model is the proper 
effective description of some real physical 2-|-l dimensional systems. Such systems should form 
an almost BPS liquid with an at least approximate realisation of the SDiff symmetry at low 
temperature. Moreover, perturbative excitations around the vacuum, governed by the relative 
strength of the potential part of the BPS baby model, on the one hand, and the non-BPS addition, 
i.e., the usual sigma model part, on the other hand, should possess rather large masses. In fact, 
there is a vast literature on planar solitons of the baby skyrme type topology in a magnetized 
medium (see, e.g., [43] and the references quoted there). The relevant topological order parameter 
is a local magnetization which can interact with an external magnetic field. It would be very 
desirable to rewrite the gauged baby BPS Skyrme model in terms of magnetization and compare 
with such systems. This issue is under current investigation. 
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